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ABSTRACT 

Properties of about 4 500 observed Ly-a absorbers are investigated using the model of 
formation and evolution of DM structure elements based on the modified Zel'dovich 
theory. This model is generally consistent with simulations of absorbers formation, 
describes reasonably well the Large Scale Structure observed in the galaxy distribution 
at small redshifts and emphasizes the generic similarity of the LSS and absorbers. 

The simple physical model of absorbers asserts that they are composed of DM 
and gaseous matter. It allows us to estimate the column density and overdensity of 
DM and gaseous components and the entropy of the gas trapped within the DM 
potential wells. The parameters of DM component are found to be consistent with 
theoretical expectations for the Gaussian initial perturbations with the WDM-like 
power spectrum. The basic physical factors responsible for the absorbers evolution are 
discussed. 

The analysis of rcdshift distribution of absorbers confirms the self consistency of 
the adopted physical model, Gaussianity of the initial perturbations and allows to 
estimate the shape of the initial power spectrum at small scales what in turn restricts 
the mass of the dominant fraction of DM particles to Mdm > 1.5 — 5 keV. Our results 
indicate a possible redshift variations of intensity of the UV background by about a 
factor of 2 - 3 at redshifts z ^ 2 - 2>. 

Key Vifords: cosmology: large-scale structure of the Universe — quasars: absorption: 
general — surveys. 



1 INTRODUCTION 

One of the most perspective methods to study the processes 
responsible for the formation and evolution of the structure 
of the Universe is the analysis of properties of absorbers 
observed in spectra of the farthest quasars. The great po- 
tential of such investigations was discussed already in Oort 
(1981, 1984) just after Sargent et al. (1980) established the 
intergalactic nature of the Ly-a forest. The available Keck 
and VLT high resolution observations of the forest provide a 
reasonable database and allow to apply statistical methods 
for such investigations. 

The essential progress achieved recently through nu- 
merous high resolution simulations of absorbers formation 
and evolution confirms that this process is closely connected 
with the initial power spectrum of perturbations. These re- 
sults allow us to consider the properties of absorbers in 
the context of nonlinear theory of gravitational instability 
(Zel'dovich 1970; Shandarin & Zel'dovich 1988) and to apply 



the statistical description of structure formation and evolu- 
tion (Demianski & Doroshkevich 1999, 2002; hereafter DD99 
& DD02) to the Ly-a forest. 

This approach is based on the modified Zel'dovich ap- 
proximation and describes the formation and evolution of 
DM structure elements for the CDM and WDM initial power 
spectra without any smoothing or filtering procedures. It 
outlines the structure evolution as a random formation and 
merging of Zel'dovich pancakes, their transverse expansion 
and/or compression and transformation into high density 
clouds and filaments. Later on, the hierarchical merging of 
pancakes, filaments and clouds forms rich galaxy walls ob- 
served at small redshifts. Main stages of this evolution are 
driven by the initial power spectrum. 

At small redshifts the main results of the statistical ap- 
proach are found to be consistent with characteristics of the 
observed and simulated Large Scale Structure (Demianski et 
al. 2000; Doroshkevich, Tucker & Allam 2002). Here we use 
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this approach for the analysis and interpretation of the ab- 
sorbers observed at large redshifts. We show that the basic 
observed characteristics of Ly-a forest are also successfully 
described in the framework of this theoretical model. 

The application of this approach to Ly-a forest was al- 
ready discussed in Demianski, Doroshkevich & Turchaninov 
(2001 a,b, hereafter Paper I & Paper II). Here we improve 
this analysis by using a richer and more refined sample of 
~ 4500 observed absorbers and a more refined model of ab- 
sorbers. Our approach uses more traditional methods of in- 
vestigation of properties of discrete absorbers rather then 
associating them with a continuous non-linear line-of-sight 
density field (see, e.g., Weinberg et al. 1998; McDonald ct al. 
2000; Croft et al. 2001; Croft et al. 2002). It allows to reach 
more clarity in the description of formation and evolution of 
absorbers and to reveal its strong link with the initial power 
spectrum. It allows also to separate several subpopulations 
of absorbers and to discuss their evolutionary history. 

The composition and spatial distribution of observed 
absorbers is complicated and at low redshifts a significant 
number of stronger Ly-a lines and metal systems is asso- 
ciated with galaxies (Bergeron et al. 1992; Lanzetta et al. 
1995; Tytlcr 1995; Lc Brunc ot al. 1996). However as was 
recently shown by Penton, Stock and ShuU (2002), even at 
small redshifts some absorbers are associated with galaxy 
filaments while others are found within galaxy voids. These 
results suggest that the population of weaker absorbers dom- 
inating at higher redshifts can be associated with weaker 
structure elements formed by the non luminous baryonic 
and DM components in extended low density regions. In 
turn, the relatively homogeneous spatial distribution of ab- 
sorbers implies a more homogeneous spatial distribution of 
both DM and baryonic components as compared with the 
observed distribution of the luminous matter. 

For the truncated Zel'dovich theory (Coles, Melot & 
Shandarin 1993; Bond&Wadslcy 1997), some statistical 
characteristics of absorbers were already discussed in Gnedin 
& Hui (1996); Hui, Gnedin & Zhang (1997); and Hui & Rut- 
ledge (1999). In particular, in these papers possible fits for 
the number density of absorbers with a threshold hydrogen 
column density, Nhi, or with a threshold Doppler parame- 
ter, 6, were proposed. However, because of the complex evo- 
lution of absorbers, these functions depend upon both the 
initial power spectrum and several random factors discussed 
below (Sees. 2.3, 4.3). So, the application and interpretation 
of these fits are limited. 

Here we use a statistical approach (DD02) based on the 
modified Zel'dovich theory which describes both linear and 
nonlinear evolution of the DM structure. For the CDM-like 
power spectrum and Gaussian initial perturbations it gives 
the expected mean ID number density and the probabil- 
ity distribution function of DM pancaJses. Both functions 
depend only on the cosmological model and moments of 
the initial power spectrum. Comparison of these functions 
with absorbers distributions confirms the self consistency of 
adopted physical model and Gaussianity of the initial den- 
sity perturbations, allows to estimate the spectral moments 
and to restrict the mass of dominant DM particles. These re- 
sults complement the investigations of the linear power spec- 
trum (see, e.g.. Croft, Weinberg, Katz & Hernquist 1998; 
Gnedin & Hamilton 2002; Zaldarriaga, Scoccimorro & Hui 
2002; Croft et al. 2002). 



However, this theory does not describe the process of 
relaxation of compressed matter, the disruption of structure 
elements due to the gravitational instability of DM pan- 
cakes and the distribution of neutral hydrogen across DM 
pancakes. Now we have only limited information about prop- 
erties of the background gas and UV radiation (Scott et al. 
2000; Schaye et al. 2000; McDonald & Miralda-Escude 2001; 
McDonald et al. 2000, 2001; Theuns et al. 2002a, b). There- 
fore, in this paper some numerical factors remain undeter- 
mined. They could be more precisely determined by special 
investigations of simulated and observed absorbers. 

This analysis should be supplemented by application of 
the discussed here approach to simulations which take simul- 
taneously into account the impact of many important factors 
and provide unified picture of the process of absorbers for- 
mation and evolution (see, e.g., Weinberg et al. 1998; Zhang 
et al. 1998; Dave et al. 1999; Theuns et al. 1999). But so far 
such simulations can be performed only in small boxes what 
restricts their representativity, introduces artificial cutoffs 
in the power spectrum and complicates the quantitative de- 
scription of structure evolution (see more detailed discussion 
in Hui & Gnedin 1997 and Paper II). The semianalitic mod- 
els (see, e.g., Bi & Davidsen 1997; Choudhury, Padmanab- 
han & Srianand 2001; Choudhury, Srianand, & Padmanab- 
han 2001; Viel et al. 2002) show other perspective approach 
to study the structure evolution at high redshifts. Perhaps, 
some results obtained below can be also useful for this pur- 
pose. 

Comparison of results obtained in Paper I and Paper 
II and in this paper demonstrates that the quality and rep- 
resentativity of the sample of observed absorbers are very 
important for the reconstruction of processes of absorbers 
formation and evolution. In this paper we use the sample of 
~ 4500 absorbers at 1.7 < z < A compiled from 18 high 
resolution spectra. However, the redshift distribution of ab- 
sorbers is inhomogcncous what increases errors and gener- 
ates some uncertainties. In particular, the statistics of ab- 
sorbers at z > 3.5 is very poor. Because of the preUminary 
character of our investigation we estimate the precision of 
only the most important parameters. 

Further progress can be achieved with richer samples 
covering the range of redshifts at least up to z --^ 4.5-5. 
The observational evidence of the reheating of the Universe 
at « w 6 (Djorgovski et al 2001; Fan et al. 2001) makes it 
important to perform a complex investigation of the early 
period of structure evolution at redshifts z ^ 4 6. 

This paper is organized as follows. The theoretical 
model of the structure evolution is discussed in Sec. 2. In Sec. 
3 the observational databases used in our analysis are pre- 
sented. The results of statistical analysis are given in Sees. 
4-6. Discussion and conclusion can be found in Sec. 7 . 



2 MODEL OF ABSORBERS FORMATION 
AND EVOLUTION 

The main observational characteristics of absorption lines 
are the redshift, Zabs, the column density of neutral hydro- 
gen, Nhi, and the Doppler parameter, 6, while the theo- 
retical description of structure formation and evolution is 
dealing with the mean linear number density of absorbers, 
nabs{z), temperature, overdensity and entropy of DM and 
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gaseous components, and with the ionization degree of hy- 
drogen. To connect these theoretical and observed parame- 
ters a physical model of absorbers formation and evolution 
is required. 

Many such models wore proposed and discussed during 
the last twenty years (see references in Ranch 1998, and in 
Paper I and Paper II). However, to describe the properties 
of the new wider and refined sample of observed absorbers 
it is necessary to develop a more detailed physical model of 
absorbers as well. Our model includes some ideas discussed 
already in earlier publications. 

2.1 Physical model of absorbers. 

In this paper we assume that: 

(i) The DM distribution forms an interconnected struc- 
ture of sheets (Zel'dovich pancates) and filaments, their 
main parameters are approximately described by the 
Zel'dovich theory of gravitational instability applied to the 
CDM or WDM initial power spectrum (DD99; DD02). The 
majority of DM pancakes are partly relaxed, long-lived, and 
their properties vary due to the successive merging and ex- 
pansion and/or compression in the transverse directions. 

(ii) Gas is trapped in the gravitational potential wells 
formed by the DM distribution. The gas temperature and 
the observed Doppler parameter, b, trace the depth of the 
DM potential wells. 

(iii) For a given temperature, the gas density within the 
wells is determined by the gas entropy created during the 
previous evolution. The g£is entropy is changing, mainly, due 
to shock heating in the course of merging of pancakes, bulk 
heating produced by the UV background and local sources 
and due to radiative cooling. 

(iv) The gas is ionized by the UV background and for the 
majority of absorbers ionization equilibrium is assumed. 

(v) The observed properties of absorbers are changing be- 
cause of merging, transversal compression and/or expansion 
and disruption of DM pancakes. The bulk heating and ra- 
diative cooling leads to slow drift of the entropy and density 
of the trapped gas. Random variations of the intensity and 
spectrum of the UV background enhance random scatter of 
the observed absorbers properties. 

(vi) In this simple model we identify the velocity disper- 
sion of DM component compressed within pancakes with 
the temperature of hydrogen and the Doppler parameter b 
of absorbers. We consider the possible macroscopic motions 
within pancakes as subsonic and assume that they cannot 
essentially distort the measured Doppler parameter. 

As compared with the model discussed in Paper II, here we 
take into account more accurately the evolution of pancakes 
after their formation. 

The formation of DM pancakes as an inevitable first 
step of evolution of small perturbations was firmly estab- 
lished both by theoretical considerations (Zel'dovich 1970; 
Shandarin & Zel'dovich 1989) and numerical simulations 
(Shandarin et al. 1995). Both theoretical analysis and simu- 
lations show the successive transformation of sheet-like ele- 
ments into filamentary-like ones and, at the same time, the 
merging of both sheet-like and filamentary elements into 
richer walls. Such continuous transformation of structure 
goes on all the time. These processes imply the existence 



of a complicated time-dependent internal structure of high 
density elements and, in particular, the unavoidable arbi- 
trariness in discrimination of such elements into filaments 
and sheets. 

Our approximate consideration cannot be applied to 
objects for which the gravitational potential and the gas 
temperature along the line of sight strongly depends on the 
matter distribution across this line. In this paper we use the 
term 'pancake' to denote structure elements with relatively 
small gradient of properties along the line of sight. With 
such a criterion, the anisotropic halo of filaments and clouds 
can be also considered as 'pancake- like'. 

The subpopulation of weaker absorbers also contains 
"artificial" caustics (McGill 1990) and absorbers identified 
with slowly expanding undcrdeusc regions (Bi & Davidsen 
1997; Zhang et al. 1998; Dave et al. 1999). Such absorbers 
produce a short lived noise, which is stronger at higher red- 
shifts z > 3. Fortunately, fraction of such absorbers does 
not exceed 10-15% of the observed sample and they cannot 
significantly distort our final estimates. 

In this paper we consider the spatially fiat ACDM model 
of the Universe with the Hubble parameter and mean density 
given by: 

H^{z) = H;',a,^{i + zf[i + nA/fim{i + zy^], (1) 

(Pm(«)) = ^a«(l + ^)^ Ho = 100/ikm/s/Mpc . 

Here Qm = 0.3 & Qa = 0.7 are dimensionless matter density 
and the cosmological term, and h = 0.65 is the dimensionless 
Hubble constant. 

2.2 Properties of the homogeneously distributed 
hydrogen 

Properties of the compressed gas can be suitably related 
to the parameters of homogeneously distributed gas, which 
were discussed in many papers (sec, e.g., Ikeuchi & Ostrikcr 
1986; Hui & Gnedin 1997; Scott et al. 2000; McDonald et 
al. 2001; Theuns et al. 2002a, b). Thus, the baryonic density 
and temperature can be taken as 

{nb{z)} = 2.4 • 10"'^cm"^(l + zf {nth^ / 0.02) , (2) 



Tig f« 1.6 • lO^KeUz), bbg = \ « 166l6gkm/s. 

V 'f^H 

Here fit, is the dimensionless mean density of baryons, 
Tftg & bbg are the temperature and Doppler parameter of the 
gas, Oig{z) ^ 1 + z describes a slow decrease of temperature 
with redshift for low density cosmological models, fcs k.mH 
are the Boltzmann's constant and the mass of the hydrogen 
atom. 

Under the assumption of ionization equilibrium of the 

gas, 

Oirnl^nHT^, ar{T)^A-10-^^(^^^y (3) 

where otriT) is the recombination coefficient (Black, 1981) 
and T~f characterizes the rate of ionization by the UV back- 
ground, the fraction of neutral hydrogen is 

Xbg = nnjub = {nb{z)) ar{Tbg)T~^ = xo(l -|- zj^, (4) 
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6.7 • 10"* fit/i^ / 16km/s 



3/2 



^° - -1^0:02 I, ' ^^ = 1°"^^^"^^- 

For the standard power law spectrum of radiation we get 

J{i') = J21 • 10~^^ (~) s~^cm~^sr~^Hz~^ , 

ri2 = 12.6J2i(3 + a^)-' . 

The gas entropy can be characterized by the function 



Tbg ^ 36keV ■ cm' / 0.02 / big{z) 



(5) 



2.3 Parameters of absorbers 

In this section wo introduce the main relations between the- 
oretical and observed characteristics of pancakes based on 
the physical model discussed in Sec. 2.1. The most impor- 
tant characteristics of DM pancakes are presented (without 
proofs) as a basis for further analysis. For more details see 
DD99 and DD02. 



2. 3. 1 Characteristics of DM component 

The fundamental characteristic of DM pancakes is the di- 
mensional, /^, or the dimensionless, g, Lagrangian thickness 
(the dimensionless DM column density) : 



(1 + .) + 

« ^h-'M^c = 33.8/z-^Mpc^^ , 



(6) 



where Iv is the coherent length of initial velocity field (DD99; 
DD02). The Lagrangian thickness of a pancake, is de- 
fined as the unperturbed distance at redshift 2 = between 
DM particles bounding the pancake. The actual thickness of 
a pancake is 



lvq{l + z) ^5 ^, 5 = Pm/{Pm{z)) ■ 



(7) 



Here 6 is the mean overdensity of a pancake above the back- 
ground density. 

Comparing the pancake surface density, n{z), with the 
distance between neighboring absorbers. 



=5.5-10^ 



Az 



H{z) 



(1 + Z)V2 



/ 0.3 , 

'—h 'Mpc, 

i In 



(8) 



we can also estimate the fraction of matter accumulated by 
absorbers as 



fab 



At(l + z) 



6.2 . 10-^gii±£l- 
^ Az 



,3/2 



0.13 



(9) 



The fraction of clustered matter characterized by {fabs{z)} 
can be compared with theoretical expectations (see Sec. 4.4). 

For Gaussian initial perturbations, the expected prob- 
ability distribution function for the DM column density is 



Ngdq ! 



1 1 
<^>"2 + ^ 



S erf(VC) 
0.82, 



3/4 



(DD02), where the parameter qo <^ 1 characterizes the co- 
herent length of the initial density field (see Sec. 7.5) and 
the dimensionless 'time' t{z) describes the evolution of per- 
turbations in the Zel'dovich theory. For ACDM model (1) 
and ^; > 2 we have 



t{z) w To 



'1 + i.2n„ 



1/3 



1 



1.27ro 



(11) 



2.2nm ) 1 + z 1 + z 

and for the amplitude of perturbations (Doroshkevich, 
Tucker & AUam 2002) 



To w (0.27 ± 0.04) 



0.2 ' 



C « 1.06g(l -h a) 



0.2 



(12) 



Strictly speaking, relations (6) and (10) are valid for 
pancakes formed and observed at the same redshift Zobs = 
Zf and after pancake formation the transverse expansion 
and compression changes its DM column density. However, 
as was shown in DD02, these processes do not change the 
statistical characteristics of pancakes observed at redshift 
^obs ^ Zf. This means that statistically we can consider each 
pancake as created at the observed redshift. More details are 
given in DD99 and DD02. 

For the DM dominated cosmological model the observed 
Doppler parameter is also closely linked with properties of 
the DM pancakes. Thus, in Paper H the 6-parameter was 
identified with the infall velocity of matter into pancakes: 



b w Vinf = uo{z)^/^'+2^, 
(z) w 2.3l^HoT^VTT^ R 



(13) 



310fem/s / 0.13 



{l + z)3/2 V ^mh"" 

This assumption is valid also some time after formation of 

the pancake but later on the pancake relaxation, compres- 
sion and/or expansion in transverse directions change the 
6-parameter. 

As is seen from (13), for (^) independent from redshift 
(10), a systematic variation of the mean Doppler parameter 
with redshift could be expected. The problem was left open 
in Paper H (due to limited observational data set) but now, 
with the more representative sample of absorbers, we see 
surprisingly weak redshift variations of the mean observed 
Doppler parameter (Sec. 4.1). This means that only small 
fraction of observed absorbers is 'young' and is described by 
(13), whereas 'older' relaxed and gravitationally confined 
absorbers dominate this sample. The simple model of re- 
laxed absorbers was already considered in Paper H. Here we 
substantially improve it. 

As is well known, for an equilibrium slab of DM the 
depth of potential well is 

.2 



(14) 



(p{z))d 

where random factors 0* & Qq characterizes the nonhomo- 
geneity of DM distribution across the slab and the evapora- 
tion of matter in the course of its relaxation. 

Analysis of numerical simulations (Demianski et al. 
2000) indicates that the relaxed distribution of DM com- 
ponent can be approximately described by the politropic 
equation of state with the power index 7 « 2 . In this case, 
we can expect that 

27 r(1.5 - 1/7) _ 4 

V^(7 - 1) r(i - 1/7) ~ 77 • 



es = 
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where r(j/) is Euler function. The actual distribution of DM 
component across a slab and the value of Q<s> depend upon 
the relaxation process which is essentially accelerated due 
to the pancake disruption into the system of high density 
clouds and filaments. This means that the factor ran- 
domly vary from absorber to absorber. 

At each pancake merging ~ 10 - 15% of matter is evap- 
orated due to the violent relaxation. This means that @q ~ 
1 for absorbers formed from the homogeneously distributed 
matter and Qq decreases for rich absorbers formed due to 
several successive mergings. 

The Doppler parameter is defined by the depth of po- 
tential well and for the iscntropic gas with pgas oc Pgil 
trapped within the well, wc get 



5 5 {p(z))S 



<5o66gy(l + 2), 



(15) 



„ 3 / Hoh 



0.13 \ / 16fcm/s 

V bbg(z) 



Variations of the gas entropy across an absorber increase the 
random variations of and 6. 



2.3.2 Characteristics of gaseous component 

The observed column density of neutral hydrogen can be 
written as an integral over a pancake along the line of sight 



Nhi = / dx ptXH = 2{xh} 



{nb{z))lyq 0.5 

1 + 2; cosd 



(16) 



Here {xh) is the mean fraction of neutral hydrogen and cos9 
takes into account the random orientation of absorbers and 
the line of sight ({cosff) ~ 0.5). Wo assume also that both 
DM and gaseous components are compressed together and, 
so, the column density of baryons and DM component are 
proportional to each other. 

However, the overdensity of the baryonic component, 
5b = nb/ (rib) is not identical to the overdensity of DM com- 
ponent, (5, (see, e.g., discussion in Matarrese & Mohayaee 
2002). Indeed, the gas temperature and the Doppler param- 
eter arc mainly determined by the characteristics of DM 
component (15) but the gas overdensity is smaller than that 
of DM component due to larger entropy of the gas. More- 
over, the bulk heating and cooling change the density and 
entropy of the gas trapped within the DM potential well. 
Under the assumption of ionization equilibrium of the gas 
(3) this process is described by the equation 



1 drib 
Ub dt 



= nbar{T)[e{T)-T-,/T\, T = Tbgh'/h, 



bg ; 



(17) 



where the recombination coefficient Qr(T) was given in (4), 
T-y ~ (5 — 10) • IO^'K characterizes the energy injected at a 
photoionization and 



e(T) w 2T4''''[1 - 0.31nT4 -l-O.lSr^/''], T4. = T/W^K, 

describes the radiative cooling for the bremsstrahlung omis- 
sion and recombination of hydrogen and helium (Black 
1981). 

The energy injected at the photoionization, T-y , depends 



1/3, 



upon the spectrum of local UV background and it varies ran- 
domly with time and space. The gas temperature varies be- 
cause of the pancake expansion and compression and due to 
the pancake disruption into a system of high density clouds. 

As is seen from (17), these processes change the baryonic 
density of pancakes and we can write 



Kb{z)5 . 



(18) 



This means that the factor K(, is small for absorbers formed 
due to adiabatic and weak shock compression because of 
the large difference between entropies of the background DM 
and the gas, and Kj, — » 1 for richer hot absorbers formed due 
to strong shock compression when this difference becomes 
small. 

Under the assumption of ionization equilibrium of the 
gas (3) and neglecting a possible contribution of macroscopic 
motions to the ^-parameter (T oc b^), for the fraction of 
neutral hydrogen and its column density we get: 

{xh) = xoKb{z)6p-^/^{l + zfe^, 13 = b/bbg , 

Nhi = Noq5/3-^^''{l + zf, iVo = 5 • lO'^m-^GH , (19) 



_ K^e^ 0.13 (cose) f Qfefe^ y / l&km/s 
~"~ ri2 Clmh^ cose I 0.02 J I bbg{z) 



3/2 



where ri2, bbg and xo were defined in (2) and (4) and the 
factor &x describes the nonhomogeneous distribution of ion- 
ized hydrogen along the line of sight. 



2.3.3 Absorbers characteristics 

Eqs. (15) and (19) link the observed and other physical char- 
acteristics of absorbers. For the DM column density, q, the 
average gas entropy, E = ln(F3), and overdensity, 5, we get: 



Nhi f/'' 
NoSo (1 + 2)« ' 



5 : 



So^{l + z) 



(20) 



exp(E) = = P'/Sl/' = K-'/'f3'/6'^' . 

Here we use the standard equation of state T/Tbg = 0^{z) = 

Fs{z)5^^'\z) ■ 

The precision of these estimates is moderate and the 
main uncertainties are generated by the unknown cosd and 
parameters 0*, Qq, @x, Kb and r-,, which vary - randomly 
and systematically - from absorber to absorber (estimates of 
q and S are independent from bbg). These variations distort 
parameters defined in (20) as follows: 



3 

q oc 



QhOi' 



.3 

^H 



Fa OC K, 



-2/3 



2/9 



(21) 



Independent estimates of these uncertainties can be ob- 
tained from the analysis of the redshift distribution of ab- 
sorbers (Sec. 6). 

However, comparing the average DM column density 
of pancakes, (q), with expectations (10) we can restrict the 
possible redshift variations of (OhBs) and estimate the com- 
bined uncertainty introduced by unknown factors. As is seen 
from (20), 



il + zfq' 



( )( 
V2- 1016 cm-2; \^ 



16km/s 



7/2 
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cos9 ( Qmh^ ^bh 



KbOxO'iQq 0.5 V 0.13 0.02 



(22) 



This relation shows that, for statistically homogeneous sam- 
ple of absorbers with {^^) — const., {cos6) = 0.5, we can 
estimate the rcdshift variations of {Gi2(z)) as follows: 

Precision of these estimates is limited because of strong ran- 
dom scatter of the product Nnih'^^ . 



2.3.4 Regular and random variations of absorber 
characteristics 

The most fundamental characteristic of absorbers is their 
DM column density, ^ « q{l + zf' . It describes the forma- 
tion and merging of pancakes, is only weakly sensitive to the 
action of random factors and defines the regular variations 
of the Doppler parameter, 6, overdensity, 5, and entropy, Fa. 
These parameters include also a random component which 
integrates the evolutionaxy history of each pancake and the 
action of random factors discussed in the previous subsec- 
tion. If the structure of a relaxed pancake can be described 
by the politropic equation of state with the effective power 
index 7 then we can discriminate the regular and random 
variations of b, 5, & Fg and introduce the reduced character- 
istics of relaxed absorbers, u, A&S: 

(l/-T-l)l 



V = ln[/?C 
A = ln[(l + z fS/f/''] = ln(5o) - 2v , 



(23) 



S' = ln[(l + z)-'F,e 



2{5/3-7)/71 



const. + lOw/3 . 



These relations indicate that, in fact, our approximate de- 
scription use only one random characteristic, namely, v, 
which is expressed through observed parameters as follows: 



■In /3- 



7- 



■In 



/ Nhi 



(24) 



67 '' 37 

As is shown in Sees. 4.5.2, the choice of a suitable value 

of 7 allows to minimize the correlation between ^ and u. Of 
course, this reduction is approximate and it can be achieved 
only statistically for a given sample of absorbers. 



2.4 Mean number density of absorbers 

Following Paper I we will characterize the ID mean number 
density of absorbers by the dimensionless function 



Ho 



H{z) dN{z) 



Ho 



dz 



(25) 



Here dN{z) is the mean number of absorbers between z and 
z + dz, {l{z)} is the mean free path between absorbers at 
redshift z and c is the speed of light. When absorbers axe 
identified with the Zel'dovich pancakes, this ID mean num- 
ber density can be linked with the fundamental character- 
istics of the cosmological model and moments of the initial 
power spectrum (DD02 & Sec. 7.5). 

As was shown in DD02, for Gaussian initial perturba- 
tions and for richer DM pancakes with sizes significantly 



larger than the coherent length of initial density field, this 
number density can be approximated by the expression 

C Wp{^thr){l + zf 



Ttabs 



HqIv 



(<j(Cthr)) 



(26) 



1 + 



nexp{^thr)[l 



WpiUr) « 0.5[1 - erf(V^)] , 

where ^thr{z) = qthr/8T'^{z), qthr{z) is the minimal (thresh- 
old) DM column density of pancakes, Wp(^thr) and {q{^thr)) 
arc the fraction of matter and the mean DM column density 
for pancakes with q > qthr and the factor (1 + «)^ describes 
the cosmological expansion of pancakes population. 

The expression (26) contains only one fitting parame- 
ter, (.thr, and it describes quite well the evolution of richer 
pancakes. However, application of this relation to evolution 
of weaker pancakes is problematic as for ^thr <C 1 we get 
from (26): 



55(1 + z) 



and this relation does not contain any fitting parameters. 

The formation of pancakes with small q is suppressed 
due to the small scale cutoff in the initial power spectrum. 
For pancakes with small qthr the redshift variations of the 
mean number density is described by the expression 



riabs 



c Vsjl + zf erf(^/e^) 
Holv 167rr(z)ygo 



exp(-6hr) : 



(27) 



(DD02) where again ^thr{z) = qthr{z)/8T^ and the value 50 
characterizes the coherent length of the initial density field 
(see Sec. 7.5). 

Both expressions, (26) and (27), are derived for the 
Gaussian initial perturbations. They are valid for DM pan- 
cakes formed due to both linear and nonlinear compressions 
and allow for merging of pancakes. Moreover, as was shown 
in DD02, if the observed redshift is identified with the red- 
shift of pancake formation then both relations remain the 
same even when the transverse compression and expansion 
of pancakes are also taken into account. This means that 
they partly account for the pancake disruption as well. 

The relation (27) describes also the impact of the 
gaseous pressure on the formation of observed absorbers. 
In this case, the parameter go must be calculated with the 
power spectrum corrected for the Jeans damping of small 
scale perturbations in the baryonic component (see Sec. 7.5). 

2.5 Observational restrictions 

The completeness of the observed samples of absorbers is re- 
stricted by the condition Nhi > Nthr ~ W^^cm~^ what in 
turn distorts the characteristics of observed absorbers and 
makes it difficult to compare them with the theoretical ex- 
pectations discussed in Sees. 2.3 and 2.4 . This means that 
these expectations should be corrected for the impact of the 
threshold column density of observed absorbers. 

Thus the investigation of spatial distribution of galaxies 
in the SDSS EDR (Doroshkevich, Tucker & AUam 2002) 
results in estimates of typical parameters of galaxy walls at 
small redshifts as 



(g)«0.4, (/3>«20, (5)«3. 



(28) 
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With these data the expected column density of neutral hy- 
drogen within the typical wall (16) is 

NhI W 0.02iVo W O.lNthrQH < Nthr . 

This moans that oven so spectacular object as the 'Greet 
Wall' does not manifest itself as absorbers. Rare absorbers 
with b > 100 km/s, Nhi ~ 10"cm~^ and ^ ~ 1.5 observed 
at 2: ~ 1.5 - 4 can be associated with embryos of wall-like 
elements of the Large Scale Structure of the Universe (see 
Sec. 5 for more detailed discussion). 

For absorbers associated with 'young' pancaJses the 
Doppler parameter is given by (13) and for the richer 'young' 
pancakes with the mean DM column density we get 

^«1, g«{g)«(l + z)-^ /3«35(l + ^)-^/^ 

S^SQs, TVhjW 10-^(1 + 2)^1/" iVoG^, (29) 

and such pancakes can be seen as absorbers already at red- 
shifts z > 1. 

For more numerous poorer pancakes with q < (q) we 

have 

30^9(1 + ^)"'^^ 5 « 4^(1 + 2)^65, (30) 
Nhi « 2.5 • 10-'iVoq^/^(l + zf^^^Qs 

-(ofe) (^) ^-^^^-' 

and such pancakes become visible for z > 1.5 only. 

These rough estimates of expected characteristics of ab- 
sorbers show that: 

(i) At redshifts z < 1.5 we can observe mainly old pan- 
cakes formed at higher redshifts which kept the measurable 
column density of neutral hydrogen up to small redshifts. 

(ii) The expected number density of observed absorbers 
rapidly increases at redshifts z ~ 1.5 - 2 when even rela- 
tively poor 'young' pancakes can pass over the observational 
threshold and become seen as weak absorbers. 

This means that our analysis can be applied to ab- 
sorbers observed at 2; > 1.5 - 2 when the impact of the 
observational threshold becomes moderate. 

2.6 Variations of the UV background 

Direct estimates of the intensity of UV background through 
the proximity effect (Scott et al. 2000) give ri2 Ri 2 ± 1 at 
2 ~ 2.5 - 3 while McDonald & Miralda-Escudc (2001) got 
ri2 ~ 0.7 — 0.4 at the same z. Indirect estimates of 
(Songaila 1998) demonstrate a probable sudden drop of the 
UV intensity of about 3 times at 2 ~ 3. This effect can be 
related to the strong ionization of Hell at these redshifts 
(Jacobsen et al. 1994; Zeng, Davidsen & Kriss 1998; Theuns 
at al. 2002a, b). 

As is shown in Sec. 6.2, possible redshift variations of 
ri2 are also seen in the redshift distribution of weaker ab- 
sorbers at z m 2.5. These variations can be fitted by the 
expression 

where ao ~ 1.5 - 3 and z^ ~ 2.5 - 3 characterize the am- 
plitude and the redshift of these variations. These estimates 
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1.5 
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2.15 


1.6 


2.1 


277 
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3.26 
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2.1 


2.7 


85 


1225 + 3179 


2.20 


1.7 


2.2 


159 


1331 + 170i° 
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1. Lu et al. (1996), 2. unpublished, courtesy of Dr. Kim 3. Hu 
et al., (1995), 4. Djorgovski et al. (2001) 5. Kirkman & Tytler 
(1997), 6. Cristiani & D'Odorico (2000), 7. Giallongo et al. (1993), 
8. Rodriguez et al. (1995), 9. Khare et al. (1997), 10. Kulkarni et 
al. (1996). 

can be distorted due to the limited representativity of our 
sample and a nonhomogeneous redshift distribution of ob- 
served absorbers plotted in Fig. 1 . 

However, weak redshift variations of the functions (6), 
(lg(A''Hj)) & (G12) (Sec. 4.1) show that the influence of variar 
tions of ri2 on the properties of absorbers is partly compen- 
sated by variations of Kh, Qx, O5&O4. introduced in Sees. 
2.3.1 & 2.3.2 . As is found in Sec. 4.2, a reasonable descrip- 
tion of the mean absorbers characteristics is achieved when 
the function (G12) (22) is 

(G12) ^G„(l- 0.172), Go = 50. (32) 

Comparison of this function with estimates of ri2 (Scott 
et al. 2000; McDonald & Miralda-Escude 2001) allows to 
estimate the unknown factors as 

(..e.e.e?)^(o.i-o.o5)(^^)\ 

More accurate estimate of these factors can be obtained with 
numerical simulations. 



3 THE DATABASE. 

The present analysis is based on the 18 spectra listed in 
Table 1. The available Ly-a lines were arranged into three 
samples. The richest sample, S}!, includes 4369 absorbers 
with lO^'^ cm"^ > Nhi > lO^^cm"^ from the first 14 high 
resolution spectra. It can be partly incomplete. For compar- 
ison, we use the most reliable sample SH which includes 
2643 absorbers from all 18 spectra for lO^^cm"^ < Nhi < 
lO^^cm"^ . These lines are more easily identified and they 
are not so sensitive to outer random infiuences. The sample 
WI4 contains 2126 weaker Unes with Nhi < lO^^cm"^ from 
the first 14 QSOs. It is mainly used to characterize possible 
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Figure 1. Redshift distribution of tlie fraction of observed ab- 
sorbers for samples W^^ (red points), Sjl (green points) and Sjg A 
(blue points). 



variations of the UV background. To decrease the scatter of 
absorbers characteristics we exclude from these samples 45 
absorbers with b > 90km/s . 

As is seen from Fig 1, for all samples the redshift distri- 
bution of absorbers is nonhomogeneous and the majority of 
absorbers are concentrated at 2 < z < 3 . This means that 
some of the discussed absorbers characteristics are derived 
mainly from this range of redshifts. Absorbers at 2; > 3.5 
were identified mainly from the spectrum of QSO 0000-260 
(Lu et al. 1996) and here the line statistics is insufficient. 



4 STATISTICAL CHARACTERISTICS OF 
ABSORBERS 

In this Section, the functions q,(,,5 & Fs axe found for 
the ACDM cosmological model (1) with btg — 16km/s, 
nth^ = 0.02 and the function {Gi2{z)) given by (32) with 



Go = 50, = {e*e,) = 0.5, 



= 0.1. 



V 



(33) 



Some quantitative results of this Section depend on the com- 
pleteness and representativity of the samples and the quite 
arbitrary choice of Qg & Oh (33) through the relations (21). 
As was noted above, for each absorber these factors change 
in the course of absorbers evolution and, in fact, the ac- 
cepted values are averaged over the sample. For the sample 
~ 10-15% of the weaker absorbers could be related to 
the artificial "noise". 

For majority of QSO spectra used in our analysis the 
scatter of observed Nhi is less then 10 per cents, while the 
scatter of the measured Doppler parameters, 6, can achieve 
20 -30 per cents. Non the less, dispersions of absorbers char- 
acteristics discussed below are defined mainly by their broad 
distribution functions and by the completeness of the sam- 
ples. Because of this, we discuss in this Section uncertainties 
of only the more interesting quantitative characteristics of 
absorbers. 




Figure 2. Redshift variations of (G12) (top panel), Doppler pa- 
rameter, (6) (middle panel) and {lg{N[ji)) (bottom panel) for 
samples SJI (blue points) and (green points). Mean values 
(34) and fit (32) are plotted by straight lines. 



samples SI4 and the mean values (6) and {IgNni) sur- 
prisingly weakly vary with redshift, 

(b) = (28.6 ± 1.6)km/s, (6) = (31.7 ± 1.6)km/s, (34) 

{IgNni) = (13.1 ± 0.22), {IgNni) = (13.6 ± 0.06), 

respectively. At the same time, the mean Doppler parameter 
systematically shifts from (b) = (25.6 ± 2.9)km/s for the 
sample Wi* up to (b) = (38.4 ± 2.7)km/s for 660 absorbers 
with Nhi > lO^'^cm"^ what indicates a weak correlation 
of observed properties of absorbers. Detailed discussion of 
observed characteristics of absorbers can be found in Kim, 
Cristiani & D'Odorico (2002), Kim et al. (2002). 

The small variations of (b) with the redshift are accom- 
panied by similar small variations of dispersion ctj, ~ 0.5{fe) 
which also do not exceed ~ 8%. This fact indicates that the 
broad distribution function of 6 is weakly dependent upon 
the redshift. The same is valid for the depth of potential 
wells, {A$) (14), and illustrates a correlation between the 
DM column density, q, and the overdensity of absorbers, S 
described by (20). The nature of such a weak redshift evolu- 
tion of these statistical characteristics of absorbers remains 
a mystery. 

The function (flu) /Go is also plotted in Fig. 2 together 
with the rough fit (32). Its weak redshift dependence is a 
natural consequence of the weak redshift variations of (b) 
and {IgNHi}- Its deviations from (32) correlate with the ab- 
sorbers distribution plotted in Fig 1 . 



4.1 Redshift variations of the mean observed 
characteristics 

Redshift variations of two observed characteristics of ab- 
sorbers, (b), and (IgNHi), are plotted in Fig. 2. For both 



4.2 Redshift variations of the mean DM column 
density, overdensity and entropy of absorbers 

Other characteristics of absorbers depend upon the physical 
model. For the model discussed in Sec. 2 (Eq. 20) with the 
function {G12) (22, 32) the redshift variations of the mean 
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Figure 3. Functions (5) = (1 + z)'^ {q) (top panel), {Fa) (middle 
panel) and the overdensity (&), (bottom panel) are plotted vs. 
redshift, 2, for samples Sj| (green points), and (blue points). 
Fits (35, 36) are plotted by lines. 



DM column density, {Q, the mean entropy and overdensity, 
{Fs) and {5), are plotted in Fig. 3. 

The mean DM column density of absorbers, (q) & (5) ~ 
{q){l + z)'^, is the most stable parameter. In principle, {^) 
does not change due to the formation and merging of ab- 
sorbers and due to their transverse compression and/or ex- 
pansion (DD02). The differences between the expected and 
measured {^) characterize, in fact, the influence of disre- 
garded factors, mainly the completeness and representativ- 
ity of the samples, the pancake relaxation and the scatter of 
the function G12. 

As is seen from comparison of Fig. 1 and Fig. 3, for both 
samples Sif and Sil variations of measured (^(2)) around 
the mean values 

(0 ~ (1- ±0.09), (C) « (0.65 ± 0.09) , (35) 

are roughly correlated with the redshift distribution of ab- 
sorbers plotted in Fig. 1 what emphasizes the limited rep- 
resentativity of the sample used in the analysis. At small 
redshifts, z < 2, the influence of the observational restric- 
tions (Sec. 2.5) increase the measured value of {Q. 

For the sample Sil the measured {£,{z)) ~ 0.65 is close 
to the expected value {^(2)) ~ 0.82 what verifies the choice 
of the amplitude of {G12) in (32). However, the estimates 
of (5(z)) depend also upon the choice of the amplitude of 
initial perturbations, ro (12), known only with a moderate 
precision. For the samples Si| , {£,{z)) increases due to the 
rejection of weaker absorbers in this sample. 

For both samples, Sil & Si| , the minimal DM column 
density is found to be qmin ~ 10~^ and it only weakly de- 
pends upon the redshift. This result is especially sensitive to 
the random noise, it is based on poor statistic and is not re- 
liable because our fit of (G12) provides a reasonable descrip- 
tion for the sample as a whole only. This value is smaUer 
than the estimate {qthr) ~ 3 • 10~^ (54) obtained in Sec. 



6.3 from the analysis of redshift distribution of absorbers. 
This disagreement requires more detailed investigation and, 
perhaps, redefinition of the factors Go, Qh and Qs . 

As is seen from Fig. 3, the redshift variations of func- 
tions {5) and {Fs) characterizing the overdensity of DM com- 
ponent and entropy of gaseous component for both samples 
Sll and Sll can be roughly fitted as follows: 

{5)«7.5r'-\ {<5)«4.7^'■^ (36) 
(F,)^1.93C'-«, C=(l + 2)/4. 

These fits emphasize the general tendencies of absorbers evo- 
lution. Comparison with the background density (2) and en- 
tropy (5) shows that the mean density of absorbers weakly 
depends on the redshift and the mean entropy weakly in- 
creases with time. These results indicate that our sample is 
dominated by long lived relaxed absorbers. 

These estimates may be rescaled according to (21). 

4.3 The mean size of absorbers 

Our model of absorbers ( see Sec. 2) allows also to estimate 
roughly the real size of absorbers along the line of sight, 
Ar. Due to the strong correlation between q and 5, this size 
should be obtained by averaging Ar as given by (7) . For the 
model parameters given in (33) and for both samples, S*!! 
and 5i|, we have 

(Ar) ~ (80 - 100)/i"^kpc , (37) 

and this size increases up to ~ 150/i^^ kpc for weaker ab- 
sorbers of the sample . For all samples the PDFs of the 
sizes are similar to the Gaussian function with a dispersion 
ar ~ 0.5(Ar). 

The mean size of absorbers in the redshift space is as- 
sociated with the Doppler parameter, b, as follows; 

(Ar,) = 2{b)/H{z) ^ 130/.-^kpc (-^Y^' ■ 

As usual, it is larger than that in the real space (37). This 
difference agrees with the domination of relaxed gravitation- 
ally bound absorbers in the sample. 

The expected mean transverse size of absorbers is com- 
parable with their Lagrangian size along the line of sight 
and was roughly estimated in DD02 as 

Difference between estimates (Ar) and (Art) agrees with 
expectations of the Zel'dovich theory and indicates that such 
absorbers could be unstable with respect to the disruption 
into a system of denser less massive clouds with Art ~ Ar 
(see, e.g., Doroshkevich 1980; Vishniac 1983). 

4.4 The mean matter fraction accumulated by 
absorbers and the mean absorber separation 

The mean matter fraction accumulated by absorbers, 
{fabs{z)), and the mean absorber separation, {Dssp{z)), de- 
fined by (8&9), are plotted in Fig. 4 for the samples W14, 
sll & -S*;!- The estimates of {fabs{z)) are sensitive to small 
measured line separations, Az ~ lO"**, which are unreli- 
able. For richer absorbers of the sample Sif, the estimates 
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Figure 4. Mean fractions of matter accumulated by absorbers, 
(faba)^ (top panel) and tlie mean absorber separations, {Dsep), 
(bottom panel), are plotted vs. redshift, z, for samples W^^ (red 
points), (green points), and SJI (blue points). Fits (39) are 
plotted by solid and dashed lines. 



of fabs(z) based on the expression (9) are also unreliable 
because they neglect the transverse compression and expan- 
sion of matter and formation of high density clouds. So, the 
representativity and precision of measured {fabs{z)) are lim- 
ited. The estimates of {Daep{z)) axe less sensitive to the first 
factor but are also influenced by the second one. The ran- 
dom scatter of the measured (fabsiz)) and {Dsep{z)) could 
be caused by the limited representativity of the samples at 
small and higher redshifts. 

In spite of uncertainties these results are interesting in 
some respects. Thus, the regular growth of {fabs{z)) and 
decrease of {Dsep{z)) with redshifts for the sample Si| indi- 
cate the progressive disruption of richer absorbers and their 
transformation into system of high density clouds. Indeed, 
due to small cross-section of such clouds, they are rarely 
cross the line of sight and their formation increases the sep- 
aration of richer observed absorber and even can decrease 
the measured (fabs) for this sample. 

Both the significant fraction of matter seen as numerous 
weaker absorbers in the sample 14^14 and its weak redshift 
variation point in favor of the domination of absorbers merg- 
ing as compared with the formation of new absorbers from 
the background matter. The merging remains unchanged the 
ifabaiz)) but progrcssively increases (Dsepiz)} with time. 
These variations can be enhanced by the observational re- 
strictions, Nhi > Nthr, discussed in Sec 2.5 . 

For all samples the fraction {fabs{z)) ~ 0.8 exceeds 
the theoretical expectation (fabaiz)) ~ 0.5-0.6 for z ~ 3 
(DD02 and Sec. 7.6) and is close to the maximal fraction, 
f,nax — 0.875, which can be compressed in the Zel'dovich' 
theory. This disagreement indicates the limited applicability 
of the one dimensional expression (9) and the limited preci- 
sion of our model of absorbers. These results must be tested 



with richer samples of observed absorbers and by compari- 
son with simulations. 

As is well known, at redshifts z > 4 the absorbers be- 
gin to overlap and their separation becomes problematic. 
Analysis of the evolution of the mean absorbers separa- 
tion allows to estimate quantitatively this effect. For this 
purpose, we can compare the mean separation of absorbers 
at various redshifts, {Daep{z)}, given by (8) with the mean 
Doppler parameter which characterizes the observed thick- 
ness of absorbers. For samples SI4 and SH these parameters 
are roughly fitted by power laws: 

/ 4 s 2.7 , r o ^2.7 

(^.ep)^(Y^j ft-'Mpc, ^[y^J h-'Mpc, 

Ho{Daep) „ / 6.1 Ho{Da.p) ^, / 8 \^'^ 

(b) ~ \l + zj ' (b) ~ \l + zj ' 

respectively. 

These results verify that at z > 4 absorbers effec- 
tively overlap and the system of individual absorbers is 
transformed into continuous absorption imitating the Gunn- 
Peterson effect. 

4.5 Distribution functions of absorber parameters 

The observed probability distribution function, PDF, of the 
DM column density of absorbers can be compared with the 
expected one (10). However, we have no theoretical expres- 
sion for the PDFs of the Doppler parameter, overdensity and 
entropy of absorbers because the action of random factors 
cannot be satisfactory described. 

Indeed, the relaxation of DM pancakes depends upon 
unknown internal structure of pancakes, the adiabatic com- 
pression and/or expansion of an absorber changes its over- 
density and temperature. On the other hand, the radiative 
cooling and bulk heating lead to the drift of the gas entropy 
and overdensity but leaves unchanged the depth of poten- 
tial well formed by DM distribution. Merging of pancakes 
increases more strongly the depth of potential well and the 
gas entropy but the overdensity of the gaseous component 
increases only moderately. All the time, the temperature 
and overdensity of trapped gas are rearranged in accordance 
with the condition of hydrodynamic equilibrium across the 
pancake. 

Because of this, our discussion of absorbers evolution 
has only phenomenological character. The PDF of the DM 
column density, Ng, is plotted in Fig. 5 for the sample 5*14 
together with the fits (10) and (40). The distribution func- 
tions of the reduced Doppler parameter, v, are plotted in 
Fig. 5 for samples Sif and 5*14 together with the Gaussian 
fits. 

4-5.1 Distribution functions of the DM column density 

The PDF for the DM column density of pancakes, Nq{q), 
plotted in Fig. 5 for samples 5*14 is the most interesting 
because it can be compared with the theoretically expected 
PDF (10) which, in particular, is sensitive to the coherent 
length of initial density field, qo- As is seen from (10), the 
PDF Nq{^), with ^ « q{l + z)^, does not depend on the 
redshift and, so, the joint PDF can be obtained for absorbers 
observed at all redshifts. 
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Table 2. Parameters of three subpopulations of absorbers for the 
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Figure 5. Top panel: 



Na, for 



the sample Fits (10) and (40) are plotted by green and blue 
lines. Bottom panel: PDF N^,, for the samples Sj| (red line) and 
SJI (green line). Gaussian fit is plotted by blue line. 



For ^ > 0.5(0, the PDF plotted in Fig. 5 is well fitted 
by the function (10) for go = and {^} = 0.57 what is close 
to {(,) = 0.65 (35). The deficit of observed absorbers with 
^ ^ {(,) can be related to incompleteness of the observed 
sample for small This result verifies the self consistency 
of the physical model used here and the Gaussianity of initial 
perturbations. 

However, as it is seen from (10), this deficit of weaker 
absorbers can be related also to the suppression of the forma- 
tion of poorer pancakes caused by the correlations of small 
scale initial perturbations and Jeans damping described by 
the parameter qo in (10). Indeed, the observed PDF is well 
fitted also by the function 



iV„ = 0.18 



rf(^) exp(-y) 
^ {1 + 0.53/y)- 



2/ = 1.15 



(40) 



what is consistent with rough estimate go ~ 0.05. This value 
depends on the statistics of absorbers at small ^ < (0 and 
can be considered as an upper limit of go only. 



4-5.2 Distribution functions of the reduced Doppler 
parameter 

As was noted above, the theoretical description for the PDFs 
of the reduced Doppler parameter, overdensity and entropy 
of the gas is problematic because they depend upon many 
random factors. Therefore, here we will restrict our analysis 
to the fits of observed PDFs, Nv, taking also into account 
the correlations of the measured b, S&iFs- 

To adjust the reduced characteristics of pancakes in- 
troduced in Eq. (23), we will minimize the correlation of ^ 
with V, A, &cS. For the sample Sif, the adjusted reduced 
characteristics are defined as follows: 
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fs is the fraction of absorbers in a subpopulation, is the linear 
correlation coefficient of b and logNui- 



A = ln[(l + zfS C^'^] = ln<5o - 2i; , 

S = ln[(l + 2)"V4"°-'^] = const. + 3.33 v . 



(41) 



These expressions correspond to Eq. (23) with the effec- 
tive power index of compressed DM component 7 — 1.8. 
For the sample SH we have 7 = 2.1. Both values are 
close to 7 = 2 discussed in Sec. 2.3.1 . With such defini- 
tion, we get the negligible correlations, < 0.01, between ^ 
and other reduced characteristics and a strong correlation 
~ 1 between v and A&iS. The distribution function A''^ 
plotted in Fig. 5 is well fitted by the Gauss function with 
(v) — 0.81 ± 0.04, f7„ ~ O.ll(u) what is consistent with the 
assumption of approximate equilibrium of majority of the 
observed absorbers. 

These results mean that, in fact, the joint action of all 
random factors is characterized by the one random func- 
tion V which can be directly expressed through the observed 
parameters as follows 



v^O.ASln l3-0.15ln{NHi /No/So), 7 = 1.8. 



(42) 



7 : 



ln(/?r 



More detailed analysis requires the discrimination of sub- 
populations of absorbers with different evolutionary histo- 
ries. 



4.6 Three subpopulations of absorbers 

As a first step of more detailed investigation of absorbers 
evolution, we can approximately separate three subpopula- 
tions of absorbers with different overdensities and entropies. 
Due to continuous distribution of both S &1F3 for the full 
samples and because both these values depend upon un- 
known random factors (21), the separation is quite arbitrary 
and characteristics of subpopulations depend upon the sam- 
ple used in the analysis and the parameters (33). Even so, 
the redshift variation of mean parameters of subpopulations 
allows to trace roughly their evolutionary history and to de- 
scribe the main stages of absorbers formation and evolution. 

The first subpopulation, let us call it "sh", with 5 > 1 
and Fa > 1 is formed due to the merging and shock compres- 
sion. The second one, called here "cl", includes low entropy 
absorbers with 5 > 1, Fs < 1 . The third one, called here 
"exp", contains peculiar absorbers with S < 1. 

This choice of the threshold parameters allows to sep- 
arate absorbers formed and observed at the same redshifts, 
Zf ~ Zobs- However, due to evolution of the background den- 
sity (2) and entropy (5), the functions 5 &iFs defined by (20) 
are shifted with time even when the overdensity and entropy 
of a long lived absorber are not changed. This shift leads to 
some artificial intermixture of subpopulations. 

For the sample Sil the mean parameters of these sub- 
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Figure 6. The mean Doppler parameter, (6) (top panel), the 
mean entropy, {Fs) (middle panel), and the mean overdensity 
{5) (bottom panel), are plotted vs. redshift, 2, together with fits 
(43, 44, &45) for the "sh" (red), "cl" (blue) and "exp" (green) 
subpopulations of the sample Sj| . 




2 3 4 



z 

Figure 7. The DM column density (top panel) and the relative 
number fractions of absorbers (bottom panel) for the "sh" (red 
lines), "cl" (blue lines) and "exp" (green lines) subpopulations of 
the sample are plotted vs. redshift, z . 



populations are listed in Table 2, where r^n is the linear 
correlation coefficient of b and Ig Nhi- Redshift variations 
of (b), (IgNni), (Fs), (5), (q), (5) and the relative num- 
ber fractions of absorbers in these subpopulations, {f,n), are 
plotted in Figs. 6 & 7 and given in (43, 44&45). 

The subpopulation "sh" is composed of the richest and 
hottest absorbers with 

{IgNni) ~ 13.0 ± 0.4, (6) ^ (37 ± 2.7)fcm/s , (43) 

(0-0.8 ±0.15, (F,) R^3.0C°-', {<5)~3.75/C, 

where ^ = (1 + 2)/4. As is seen from Fig. 7, for this sub- 
population the decrease of the relative number fraction of 
absorbers, (fm), with z is accompanied by the growth of 
the DM column density, (q). Such variations are typical for 
absorbers formed due to pancake merging and strong shock 
compression of the matter. 

The subpopulation "cl" includes absorbers with 

(/ffA^H/) ~ 13.5 ± 0.37, (&) ^ (23.1 ± l.l)km/s, (44) 

(0-0.69 ±0.1, (F,) ^ 0.6 ±0.08, (5) ^ 8.7C"^■^ 

Such behavior of the mean characteristics indicates that this 
subpopulation can be dominated by absorbers formed in low 
entropy regions at ~ Zots due to adiabatic and weak shock 
compression. Such compression provides higher overdensity 
than the strong shock compression and does not change the 
entropy of compressed matter. For these absorbers the ra- 
diative cooling can also be more important. This subpop- 
ulation includes also some fraction of long lived absorbers 
formed at > Zobs due to merging and strong shock com- 
pression. This intermixture is caused by redshift variations 
of background parameters (2) and (5). 

The peculiar subpopulation "exp" accumulates poorer 



absorbers with 

{IgNHi) ~ 12.3 ± 0.3, (6) ^ (18 ± 5.4)km/s, (45) 

(g)^l.l■10~^ (Fs)«2.1C^■^ ((5) ^ 0.6 ± 0.08 . 

Such behavior of (q), (S) & (Fs) shows that these absorbers 
could be related to the unstable poorer pancakes formed 
within low density regions (see, e.g., Zhang et al. 1998). 
They are observed mainly at higher redshifts, they expand 
together with the background and disappear at low redshifts. 
Large values of the factor 7 = 3.5 and of the reduced Doppler 
parameter (v) — 1.3 (Table 2), verify the peculiar character 
of this subpopulation. 

Moreover, the analysis performed in Sec. 6.3 indicates 
that, for weaker absorbers dominating this subpopulation, 
at least the DM column density, (q) , can be underestimated. 
This means that other parameters of the subpopulation can 
also be estimated unreliably. In particular, if these absorbers 
are unstable then the random factors become time depen- 
dent and the estimates (45) must be corrected. Fortunately, 
the fraction of such absorbers is small and its influence on 
other estimates is limited. 

In spite of the quite arbitrary separation of subpop- 
ulations these results illustrate the action of evolutionary 
factors mentioned above. The objective character of the dis- 
crimination is seen from the comparison of observed param- 
eters, {IqNhi) and (6), and their linear correlation coeffi- 
cient, TbH, listed in Table 2 for the full samples and for the 
separated subpopulations. 

These results confirm that majority of observed ab- 
sorbers represent gravitationally bounded pancakes formed 
in the course of adiabatic and shock compression at 2:/ > 
Zobs- They demonstrate similarity of evolutionary history 
and generic origin of subpopulations "sh" and "cl". More 
detailed statistical investigation of absorbers evolution re- 
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quires finer discrimination of subpopulations and tlierefore 
richer set of observed absorbers is needed. 



4.7 Absorbers as a test of background properties 

Some absorbers characteristics can be used to estimate the 
redshift variations of mean properties of homogeneously dis- 
tributed hydrogen (see, e.g., Hui & Gnedin 1997; Schaye ct 
al. 1999, 2000; McDonald et al. 2001). Thus, weak redshift 
variations of (Fs) (44) suggest also weak redshift variations 
of Tbg (2). However, such estimates are inevitably approx- 
imate and their significant scatter is caused by action of 
many factors discussed above. 

To obtain more stable results Schaye et al. (1999, 2000) 
consider a cutoff at small 6 in the distribution of b{NHi)- 
Indeed, relations (20) can be rewritten as follows: 

*=K,»r*:"(^)°"*(^)"'. («) 

Parameters A'o, So, bbg & Kb were defined by (2, 15, & 19). Ev- 
idently, for Fs = const, this relation is identical to the equa- 
tion of state = const. 5^''''. It demonstrates that, for the 
subpopulation of absorbers formed at a/ « Zobs with the 
same Fs, the observed parameters z, b and Nhi are strongly 
correlated and, for a given Nhi, b decreases with Fg. How- 
ever, statistics of absorbers near the cutofT is small, the ac- 
tion of factors discussed above erodes the cutoff and makes 
its interpretation less reliable. 

For absorbers selected from the sample SI4 in four red- 
shift intervals, z < 2 < z < 2.5 < z < 3 < z, the boundary 
of the observed distribution b{NHi) can be roughly fitted by 
the relation 

b^boil + z)~\NHi/W^'^cm~'^y'' . (47) 

For 2 < 3 we have 75 ~ 0.25 that coincides with (46) but the 
value of bo ~ 22 km/s is ~ 2 times smaller, what is expected 
from (46). This means that in the framework of the model 
under consideration the cutoff can be naturally related to a 
fraction of low entropy absorbers with Fs ~ 0.3 rather then 
with Fs ~ 1. These absorbers can be formed at 2/ ~ Zobs 
due to the adiabatic compression within low entropy regions 
or due to both adiabatic and shock compressions at redshifts 
Zf > 1.5 Zobs- For richer absorbers with IqNhi > 13 the 
radiative cooling is also more important. This explanation 
is quite consistent with the estimates (44). 

At 2 > 3 we have for the fit parameters 7(, ^ 0.15 
0.17, bo ~ 38km/s, what differs from (46). This difference 
can be partly related to the poor statistics at higher redshifts 
and to admixture of long lived absorbers formed at redshifts 
2 > 4 when stronger systematic and random variations of 
UV background are expected. In particular, at 2 ~ 3 Hell 
becomes ionized and the spectrum of UV background and 
the effective power-low index of the temperature-density re- 
lation are changed (Songaila 1998; Schaye et al. 2000; The- 
uns et al. 2002a, b). 

These results demonstrate that, in the framework of 
model under consideration, the cutoff in the distribution 
of b(NHi) can be explained but its interpretation is not so 
clear. The problem requires further investigation. 



5 WALLS AT HIGH REDSHIFTS 

Majority of the rare absorbers with higher Doppler param- 
eter, b > (80 - 90)km/s and moderate column density of 
neutral hydrogen Nhi < lO'^'^'^cm"^, can be identified with 
the embryos of richer structure elements which are seen now 
as rich walls in the observed galaxy distribution. The num- 
ber of such absorbers is small and their statistical char- 
acteristics cannot be reliably determined. However, for 45 
such absorbers with b > 90km/s selected in 9 spectra at 
1.8 < 2 < 3.8 we have 

{IgNni} ~ 13.2, (0 ~ 1-45, (<5) ~ 6, (F,)~10, (48) 
and the mean comoving separation of such absorbers (8) is 
(Ds.p) « (60±20)V0.3/n^/i-'Mpc. (49) 

Such large separation suggests that these objects will 
retain their individuality up to small redshifts and, indeed, 
this separation is quite consistent with the mean wall sep- 
aration (Dw) ~ (66 ± I3)h~ Mpc measured for the SDSS 
EDR at 2 = (Doroshkevich, Tucker & Allam 2002). How- 
ever, for smaller b, the number of absorbers rapidly increases 
and already for 139 absorbers with b > 70km/s in 15 spectra 
the mean separation decreases to Dsep ~ (40 ± 20)/i~^Mpc. 

These results show that the walls begin to form already 
at 2 ~ 3 and continue to form up to the present time. At 
the same time, the continuous distribution of all parameters 
of absorbers shows that at high redshifts the discrimination 
of walls and other structure elements is quite arbitrary. The 
problem deserves further investigation in a wider range of 
redshifts with a more representative sample of absorbers es- 
pecially at high redshifts. Perhaps, such walls can be also 
observed in the galaxy distribution at high redshifts. 



6 REDSHIFT DISTRIBUTION OF 
ABSORBERS 

In this section we compare the observed redshift distribution 
of absorbers with theoretical relations (26, 27) describing the 
expected redshift evolution of the ID number density of DM 
pancakes, Uabs, due to their formation and merging. As is 
seen from these relations, the number density depends upon 
the cosmological model, moments of initial power spectrum 
and the threshold DM column density, qthr- Therefore, to 
compare these expectations with ol)S('rve<l absorbers distri- 
bution we have to select samples with q > qthr rather than to 
use samples with a minimal Nhi discussed in Sec. 4. Such 
samples can be prepared using estimates of q as given by 
(20) in spite of their unreliability for smaller q (see Sec. 4.2, 
6.3). 

Both relations (26&27) are derived for the Gaussian 
initial perturbations. They depend on the redshift of ab- 
sorbers only and, so, the impact of absorbers velocities and 
Doppler parameters on the resulting estimates is minimal. 
At 2 < 4 these factors increase the scatter of measured Uabs 
but only at 2 > 4 - 4.5 they lead to a significant overlapping 
of absorbers (Sec. 4.4). The Jeans damping caused by the 
temperature of background is well described by (27) through 
the parameter go. 

Four samples, namely, Qi", Q14, QI4, and QH, with 
qtf^^ = 0.05, 0.02, 0.01 & 0.001 were selected from 14 high 
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resolution spectra. These samples contain 1554, 3299, 3998 
and 4475 absorbers, respectively. Comparison of qthr used 
for the sample preparation with estimates of (qthr) obtained 
from fits (26, 27) allows to test the self consistency of our 
physical model and the choice of parameters (33) . 

The richest sample includes almost all lines and 
probably is incomplete. It can be used for comparison with 
results obtained for other samples. The sample Qil contains 
richer absorbers and is weakly sensitive to properties of small 
scale density field described by the parameter qo in (27). It 
is used to test the expression (26) and to estimate (qthr) 
for such absorbers. The samples Q14, Q14 can be used for 
independent estimates of both (qthr) and go - Of course, these 
estimates are statistical and averaged over the sample. 

To test the sample dependence of our results we use the 
sample Qig with qthr ~ 0.01 selected from all 18 spectra. As 
compared with the sample Q14 it contains ~10% additional 
absorbers with Nhi > 10^'^ cm~'^. To demonstrate the pos- 
sible variations of intensity of UV background we consider 
also the sample Wii containing 2126 weaker absorbers with 
Nhi < lO^^cm"^. Here we also use the sample of HST data 
(Bahcall et al. 1993, 1996; Jannuzi et al. 1998) which con- 
tains 1000 absorbers at z < 1.5 . This sample is probably 
incomplete at least at 2; > 1 . 



6.1 Selection of Poissonian subsamples of 
absorbers 

The first problem is the selection of Poissonian subsamples 
of absorbers and estimation of their mean number density at 
various redshifts. For this purpose, we use the measured red- 
shift separations of neighboring absorbers, Az, what partly 
attenuates the influence of selection effects inherent in indi- 
vidual spectra. Instead of separation we use the dimension- 
less comoving distance, 

Al = HoAz/H{z), 

and absorbers in the interval z — dz < z < z + dz, taken 
from the sample under investigation, are organized into an 
'equivalent single field' by arranging the separations Al one 
after the other along the line of sight. For richer samples, 
distribution of absorbers separations obtained in this way is 
similar to the Poissonian one and the mean number density, 
(n) , can be found by comparing the measured PDF with the 
differential or cumulative Poissonian PDFs: 

dN/dx = exp(-{n) x)/{n) , (50) 

N = exp{-{n) x) . (51) 

To decrease uncertainties and to estimate the actual 
scatter of these fits we rejected all points before the maxima 
of differential PDF and only more representative middle part 
of PDFs with the fraction of points / < fthr were used. 
The threshold fraction, fthr, was varied between fthr ~ 0.7 
and fthr ~ 0.95. Mean values and dispersions of the set of 
measurements with different fthr were taken as the actual 
value and scatter of nabs. As a rule, the fit (51) of cumulative 
PDFs gives more stable estimates and smaller error bars 
than the fit (50) of differential PDFs. Both estimates of Uats 
are plotted in Figs. 8&9. 
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Figure 8. Top panel: redshift distribution of absorbers, riabsi^), 
found with fits (50) and (51) (red and green points) for the sample 
Fit (27) is plotted by green line. Bottom panel: redshift 
variations of rtabs/i^fit for the same sample. Fit (52) is plotted 
by blue line. 

6.2 Variations of intensity of the UV background 

For the sample W14 the redshift distribution of absorbers 
is plotted in Fig. 8 together with the fit (27) which charac- 
terizes the expected redshift distribution of absorbers. The 
weak absorbers in this sample are especially sensitive to vari- 
ations of intensity of the UV background radiation what is 
manifested by the strong variations of riabs in comparison 
with the smooth fitting curve. This effect can be enhanced 
by the incompleteness of the sample and the nonhomogene- 
ity of redshift distribution of observed weaker absorbers. At 
z < 2 the rapid decrease of Uats is caused by the influence 
of the threshold of the observed column density of neutral 
hydrogen, Nhi > Nthr ~ W^^cm~^, discussed in Sec. 2.5. 

The observed variations of (uabs) are well fltted by the 
function 

{nabs)/nfu = 1.6 + 0.65 th (^^-j^) , (52) 

also plotted in Fig. 8. Similar variations with smaller ampli- 
tude are also seen for absorbers in other samples. 

The column density of neutral hydrogen is proportional 
to and an increase of F.^ shifts weak absorbers below the 
observational threshold. This means that the relation (52) 
can describe quite well the redshift dependence of variations 
of the UV background but their amplitude depends upon the 
distribution function of Nhi- Direct estimates show that the 
variations of the intensity by about 2-3 times described by 
(31) are consistent with the observed variation of (uabs)- As 
was noted in Sec. 2.6 the influence of these variations on the 
mean characteristics of absorbers can be partly compensated 
by variations of random factors Qh & 9,5. 

Variations of the UV background at 2 ~ 3 were al- 
ready detected by different methods (see, e.g., Songaila 1998; 
Scott et al. 2000; McDonald & Miralda-Escude 2001). This 
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Figure 9. Top panel: redshift distribution of absorbers, nabai^), 
found witli fits (50) and (51) (red and green points) for the sample 
QJI^. Fit (27) is plotted by the blue line. Bottom panel: DilTerences 
between measured and fitted absorber distributions. 



approach seems to be sufficiently perspective but quantita- 
tive results are sample dependent and must be tested with 
more representative and homogeneous samples of weak ab- 
sorbers. 



6.3 Redshift distribution of DM pancakes 

In this Section we consider the redshift distribution of DM 
pancakes for the samples Q14, Qil, Qlt, and QH prepared 
with qthr(z) ~ const, for parameters (33). For the sample 
(5i4, the function nabs{z) is plotted in Fig. 9 together with 
the best fit (27) and the random scatter of observed points 
around the fit. At redshifts 1.8 < z < 3.3 where the repre- 
sentativity of samples is better the scatter does not exceed 
~ 20% . 

At redshifts z ~ 1.5 - 2 the rapid growth of observed 
mean number density of absorbers 



(53) 



coincides with the expected one (30) and is naturally ex- 
plained by the influence of observational threshold Nhi > 
Nthr ~ lO^^cm-^. Such cutoff is not seen for the subpopula- 
tion of richer absorbers with Nhi > lO^^cm"'^ what verifies 
its connection with the observational threshold. The proba- 
ble incompleteness of the sample of weak absorbers ai z < 
1.5 enhances this rapid growth of {uabs)- 

The best estimates of the fit parameters (27) for the 



samples Q14, 



)10 

;14, 



and Qi4 are, respectively, 



(qthr) ~ (3.3 ± 0.3) ■ 10"", (qthr) ~ (3.1 ± 0.3) • 10' 

{?^^r)~(3.1±0.3)•10-^ 
and for all three samples 
go « (0.9 ±0.1) ■ 10"^ . 



(54) 



(55) 



and the resulting estimates 



,-2 



Here the formal errors of the fits are given. 

To test the sample dependence of the results we ana- 
lyzed in the same manner the sample Sls selected with qthr — 
0.01 from all 18 spectra (4404 absorbers). In this case, the 
excess of absorbers with Nhi > 10^'^ cm"'^ and larger separa- 
tions decrease the measured riabs 
are 

(qthr) ~ (3.2 ± 0.3) X 10-^ go ~ (1-1 ± 0.12) x 10 

We see that such perturbations lead to moderate variations 
of both (qthr) and go ■ 

Real precision of both estimates (54) and (55) depends 
upon several factors, mostly on the limited representativity 
of the sample and precision of estimates to (12). As is seen 
from (27), we really evaluate the product qthrTQ^ and goT"o , 
and, so, the real precision of both estimates (54) and (55) is 
not better than ~ 20 - 30% . 

As was noted in Sec. 2.4, the redshift distribution of 
absorbers with larger qthr is weakly sensitive to small scale 
correlations of perturbations and it can be fitted by the ex- 
pression (26). Indeed, for samples Qf2 a-nd Q14 these dis- 
tributions are well fitted by expression (26) with (qthr) ~ 
0.035 & 0.05, respectively, and the amplitude of the fit is 
only ~ 1.2 - 1.3 times larger then what is expected in (26). 
However, already for the sample Q\1 the difference between 
expected and measured amplitudes of the fit (26) increases 
by up to ~ 1.75 times. 

These results show that the actual reliability and pre- 
cision of our estimates (55) are limited due to limited repre- 
sentativity of the samples and limited precision of estimates 
of the amplitude to (12). None the less, they suggest that: 

(i) The redshift distribution of absorbers is quite well fit- 
ted by relations (26) and (27) derived for the Gaussian initial 
perturbations. 

(ii) The distribution of weaker absorbers is actually sen- 
sitive to the small scale cutoff of power spectrum described 
by the parameter go ~ (0.9 ± 0.3) • 10"^ . 

(iii) Comparison of estimates (54) with qthr used for the 
sample selection indicates that, for weak absorbers with 
g ~ 10"^, the expression (20) with parameters (33) un- 
derestimates the DM column density. However, for richer 
absorbers the divergence decreases to the quite moderate 
value (~ 10%). 

The model as a whole and estimates of go can be im- 
proved with richer sample of absorbers especially at high 
redshifts z > 3.5 where the representativity of our samples 
is insufficient. However, as was noted in Songaila (1998) and 
Schaye et al. (2000), at z > 3 the random variations of the 
UV background and the scatter of (riabs) increase. 



7 SUMMARY AND DISCUSSION. 

In this paper we continue the analysis initiated in Pa- 
per I and Paper II based on the statistical description of 
Zel'dovich pancakes (DD99, DD02). This approach allows 
to connect the observed characteristics of absorbers with 
fundamental properties of the initial perturbations without 
any smoothing or filtering procedures, to reveal and to illus- 
trate the main tendencies of structure evolution. It demon- 
strates also the generic origin of absorbers and the Large 
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Scale Structure observed in the spatial galaxy distribution 
at small redshifts. 

We investigate the new more representative sample of 
~ 4 500 absorbers what allows us to improve the physical 
model of absorbers introduced in Paper I and Paper II and 
to obtain reasonable description of physical characteristics 
of absorbers. The progress achieved demonstrates again the 
key role of the representativity of observed samples for the 
construction of the physical model of absorbers and reveals 
a close connection between conclusions and observational 
database. 

However, the representativity of this sample is not suf- 
ficient for the more detailed study of structure evolution. 
Further progress can be achieved with richer sample of ob- 
served absorbers. 

7.1 Main results 

Main results of our analysis can be summarized as follows: 

(i) The approach used in this paper allows to link the ob- 
served and other physical characteristics of Ly-a absorbers 

such as the ovcrdcnsity and entropy of the gaseous compo- 
nent and the column density of DM component accumulated 
by absorbers. 

(ii) The basic observed properties of absorbers are quite 
successfully described by the statistical model of DM con- 
fined structure elements (Zel'dovich pancakes). Comparison 
of independent estimates of the DM characteristics of pan- 
cakes confirms the self consistency of the physical model for 
richer absorbers. However, some characteristics of pancakes 
formation and evolution remain uncertain. 

(iii) In the framework of this approach, all characteristics 
of absorbers can be expressed through two functions, 
describing the systematic and random variations of absorber 
properties. They are directly expressed through the observed 
parameters, z, NHiSzb. 

(iv) The main stages of structure evolution arc illustrated 
by separation of three subpopulations of absorbers with high 
and low entropy and low overdensity. 

(v) The absorbers with high Doppler parameter, 6 > 90 
km/s, can be naturally identified with the embryos of wall- 
like structure elements observed in the spatial distribution 
of galaxies at small redshifts. 

(vi) The strong suppression of the mean number density 
of absorbers at z < 1.8 can be naturally explained by the 
influence of the observational threshold Nhi > lO^^cm^'^. 

(vii) Redsliift distribution of weaker absorbers indicates 
the probable systematic redshift variations of the UV back- 
ground. They can be caused by the reionization of Hell. 

(viii) The PDF of the DM column density and the redshift 
distribution of absorbers are consistent with Gaussian initial 
perturbations. 

(ix) We estimate the spectral moment mo which in turn 
is linked with the cutoff of the power spectrum caused by 
the mass of dominant fraction of DM particles and the Jeans 
damping. 

7.2 Variations of intensity of the UV background 

The intensity of the UV background and its variations at 
redshifts 2: ~ 2 - 3 were detected by different methods (see. 



e.g., Songaila 1998; Scott et al. 2000; McDonald & Miralda- 
Escude 2001). However, the achieved precision of these esti- 
mates is limited. The variations of background temperature 
at these redshifts were also detected by Schaye et al. (2000) 
and McDonald et al. (2001). 

The redshift variations of the mean number density of 
weak absorbers discussed in Sec. 6.2 can be naturally related 
to such variations caused by the reionization of Hell and pos- 
sible variations of activity of quasars and other sources of 
UV radiation. They can be enhanced by the incompleteness 
and limited representativity of our samples and by the non- 
homogeneous redshift distribution of observed absorbers. 

However, these variations do not lead to appreciable 
variations of the mean absorbers characteristics discussed 
in Sec. 4 what indicates the complex character of absorbers 
evolution. The problem requires more detailed investigation 
using both observations and simulations. 

7.3 Properties of absorbers 

The physical model of absorbers introduced in Sec. 2 links 
the measured z, b and IqNhi with other physical charac- 
teristics of both gaseous and DM components forming the 
observed absorbers. Uncertainties in the available estimates 
of the background temperature and UV radiation and un- 
known parameters of the model restrict its applications. For- 
tunately, actions of these factors partly compensate each 
other, what allows us to obtain reasonable statistical de- 
scription for majority of absorbers. The self consistency of 
this approach is confirmed by two independent estimates 
of the column density of DM pancakes. Numerical simula- 
tions can clarify action of unknown factors and improve the 
model. 

Analysis of mean absorbers characteristics performed 
in Sec. 4 shows that the sample of observed absorbers is 
composed of pancakes with various evolutionary histories. 
We discuss five main factors that determine absorbers evo- 
lution after formation. They are: the transverse expansion 
and compression of pancakes, their merging, the radiative 
heating and cooling of compressed gas, and the disruption 
of structure elements into a system of high density clouds. 
The first two factors change the overdensity of DM and gas 
but do not change the gas entropy. Next two factors change 
both the gas entropy and overdensity but do not change the 
DM characteristics. Using the entropy and overdensity we 
can roughly discriminate three subpopulations of absorbers, 
illustrate the influence of these factors and correlations be- 
tween absorbers characteristics. 

Introduction of the reduced Doppler parameter, v, over- 
density, A, and entropy, S, allows to discriminate between 
the systematic and random variations of absorbers prop- 
erties. The former ones are naturally related to the pro- 
gressive growth with time of the DM column density of ab- 
sorbers, q[z), what can be described theoretically. On the 
other hand the action of random factors cannot be sat- 
isfactory described by any theoretical model. However, in 
the framework of our approach, the reduced parameters are 
found to be strongly correlated (41) and, in fact, the joint 
action of all random factors is summarized by one random 
function, v, directly expressed through the observed param- 
eters (42). These results alleviate the problem of absorbers 
description and, perhaps, the modeling of the Ly-a forest 
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based on the simulated DM distribution (Viel et al. 2002). ^ 
However, estimates of the matter fraction accumulated by^' 
absorbers (Sec. 4.4) demonstrate that the model needs to be 
improved. ^ 

Estimates of the size of absorbers (Sec. 4.3) suggest a 
possible fast disruption of richer absorbers into a system of 
high density clouds, what is the first step of formation of 
dwarf galaxies. ^ 

o 

7.4 Absorbers as elements of the Large Scale 
Structure of the Universe 

The close connection of absorbers with the Large Scale 
Structure (LSS) of the Universe observed in the spatial ^ 
galaxy distribution was demonstrated in numerical simu-^-^ 
lations and was recently confirmed by direct observations ^ 
(Penton, Stocke & ShuU 2002). Our results indicate the 
generic link of absorbers and DM Zel'dovich pancakes and 
demonstrate that the absorbers characteristics are consis- 
tent with Gaussian initial perturbations and the Harrison- 
Zel'dovich power spectrum. The possible identification of 
some absorbers with embryos of walls observed in the galaxy 
surveys at small redshifts shows that, perhaps, also a wider 
set of richer absorbers could be identified with the LSS ele- 
ments at high redshifts. 



7.5 Characteristics of the initial power spectrum 

The amplitude and the shape of large scale initial power 
spectrum are approximately established by investigations of 
relic radiation and the structure of the Universe at 2 ^ 1 
detected in large redshift surveys such as the SDSS (Dodel- 
son et al. 2002) and 2dF (Efstathiou et al. 2001). The shape 
of small scale initial power spectrum can be tested at high 
redshifts where it is not so strongly distorted by nonlinear 
evolution (see, e.g.. Croft et al. 2002). 

Recent discussions on the mass of dominant fraction of 
DM component and the shape of small scale initial power 
spectrum are focused on the formation of low mass halos 
in CDM simulations in comparison with observed low mass 
satellites (see, e.g., Colin, Avila-Reese & Valenzuela 2001; 
Bode, Ostriker & Turok 2001), and at the simulations of 
absorbers formation. Thus, Narayanan et al. (2000) esti- 
mate the low limit of this mass as Mdm > 0.75 keV, while 
Barkana, Haiman & Ostriker (2001) increase this limit up 
to Mdm > 1 - 1.25 keV. 

Here we can compare these estimates with direct mea- 
surements of redshift distribution of absorbers. In fact, our 
estimates are related to two spectral moments of the ini- 
tial power spectrum, m_2 and mo. The first of them de- 
pends mainly upon the large scale power spectrum and was 
independently estimated from the measured characteristics 
of observed galaxy distribution at small redshifts. Together 
with cosmological parameters, Q.m and h, this moment de- 
fines the coherent length of initial velocity field, (6), 
used in Sec. 2.3 for discussion of properties of DM distri- 
bution. The second moment, mo, depends upon the small 
scale power spectrum and the shape of transfer function. 
This means that the estimates of the mass of DM particles 
are model dependent even though our estimates of go and 
spectral moments are based on the analysis of observations. 
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Figure 10. Spectral moments, m_2 and mo, and the parame- 
ter go E're plotted vs. Mj^m for the CDM (red points) and WDM 
(blue points) models with the damping factors Dw (57) and with 
both damping factors, Dw^Dj (58) for xj = (2&0.5) • IQ-^, 
(green points). Black lines show the observational restrictions 
(59). 



As was shown in DD99 and DD02, for the CDM and 
WDM models with the Harrison - Zel'dovich asymptotic 
of power spectrum, p{k) oc k, and transfer functions, T{k), 
given in Bardeen et al. (1986) the coherent length of ini- 
tial density perturbations, Ip, and the parameter go can be 
written as follows: 



go = 5 — —, In = qolv ~ 0.34^~^Mpc ( 
mo V 



go 0.2 
0.01 ^rah 



(56) 



m-2 = / dx xT'^{x)^QmZ, mo = dxx^T^{x), 
10 Jo 



Mpc, 



X = fco = Qmh^ /Mpc, /„ = ^ = ^^'^„ 

ko ko i/m_2 ilm.h^ 

where k is the comoving wave number. 

For WDM particles the dimensionless damping scale, 
Rf, and the damping factor, Dw, are 

\ Mdm J 

Dw = exp[-xRf - [xRff] , (57) 

(Bardeen et al. 1986). The Jeans wave number, kj, and the 
damping factor, Dj, can be taken as 

kj"- ^0.7bbg{l + z)H-\z), Dj ^ {1 + x^x'')-^ (58) 

4 bbg 



xj = ko/kj « 2 ■ 10" 



0.13 \l 1 + z I6km/s 



(see, e.g., Matarrese & Mohayaee 2002). 

Variations of the spectral moments and go with the 
mass of DM particle obtained by integration of the power 
spectrum with these damping factors are plotted in Fig. 
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10 for the CDM and WDM transfer functions and for 
= 2 ■ 10"^ & 0.5 • 10"^. As is seen from this Fig., the 
low limit go > 0.07^/a^ caused by the Jeans damping (58), 
restricts the range of masses of DM particles detectable with 
this approach to Mdm < (3 - 5)keV. 
Our la estimates 

go ~ (0.6 - 1.2) ■ 10"^ mo -0.15 -0.5, 

Mdm « (1.5- 5)keV, (59) 

are close to those of Narayanan et al. (2000) and Barkana, 
Haiman & Ostriker (2001). They are based on absorbers 
with the redshifts z < 3.5 where the scatter of measured 
Uabs is minimal and, so, they weakly depend upon both the 
Doppler parameter and velocities of absorbers which essen- 
tially distort Uabs at higher redshifts (Sec. 4.4). 

There are plenty of possible candidates for the WDM 
particles with the mass ~ 1 keV. They are, for example, the 
sterile neutrinos, majorons and even shadow particles. De- 
tailed review of possible candidates can be found in Sommer- 
Larsen & Dolgov (2001) and Dolgov (2002). However, our 
estimates use the spectral moments and therefore do not 
forbid the existence of multicomponent heavy DM particles 
with both thermal and non-thermal (Lin et al. 2001) distri- 
butions with the same spectral moments m_2 and mo. 

7.6 Reheating of the Universe 

Recent observations of high redshift quasars with z > 5 
(Djorgovski et al. 2001; Becker et al. 2001; Pentericci et al. 
2001; Fan et al. 2001) demonstrate clear evidences in fa- 
vor of the reionization of the Universe at redshifts 2; — 6 
when the volume averaged fraction of neutral hydrogen 
is found to be fu > 10^^ and the photoionization rate 
r.^ ~ (0.2 - 0.8) • lO^^^s"^ . These results are consistent 
with those expected at the end of the reionization epoch 
which probably takes place at z ~ 6. Extrapolation of the 
mean separation of absorbers discussed in Sec. 4.4 is consis- 
tent with these conclusions. 

These results can be compared with expectations of the 
Zel'dovich approximation (DD02) for go = 10~^. The poten- 
tial of this approach is limited since it cannot describe the 
nonlinear stages of structure formation and, so, it cannot 
substitute the high resolution numerical simulations. How- 
ever, it describes quite well many observed and simulated 
statistical characteristics of the structure such as the red- 
shift distribution of absorbers and evolution of their DM 
column density. This approach does not depend on the box 
size, number of points and other limitations of numerical 
simulations and it successfully augments them. 

This approach allows to estimate the fractions of DM 
component accumulated by high density clouds, fc.i, fila- 
ments, //, and pancakes, /pan, at different redshifts. These 
functions plotted in Fig. 11 for go = lO"'^ show that at z ~ 6 
only ~ 3.5% of the matter is condensed within the high den- 
sity clouds which can be associated with luminous objects. 
This value can increase up to ~ 5 - 6% with more correct 
description of the clouds collapse. At the same redshifts, 
~ 27% and ~ 12% of the matter can be already accumulated 
by pancakes and filaments, respectively. These expectations 
are smaller than estimates of (fabs) obtained in Sec. 4.4. 

The Zel'dovich approximation allows also to estimate 
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Figure 11. Top panel: expected redshift variations of the DM 
fraction accumulated by high density clouds, filaments and pan- 
cakes (red, green and blue lines), respectively, at go = 10~^. Bot- 
tom panel: the expected mass function of the DM clouds at z = 6 
and go = 10~^. 

the mass function of all structure elements (DD02) at dif- 
ferent redshifts. For go — 10~^ and at z ~ 6, this function is 
plotted in Fig. 11. For these go and 2, the mean DM mass of 
structure elements is expected to be ~ 10^^ Mq and the main 
mass is concentrated within clouds with Md ~ 0.1{Afc;). As 
is seen from Fig. 11, majority of the clouds have mass be- 
tween 10~'^{Mci) and 10 (Md). The formation of low mass 
clouds with Mci < 10^ M© is suppressed due to strong cor- 
relation of the initial density and velocity fields at scales 
< Zp ~ 0.3/i"^ Mpc (56). However, the numerous low mass 
satellites of large central galaxies can be formed in the course 
of disruption of massive collapsed clouds at the stage of their 
compression into thin pancake-like objects (Doroshkevich 
1980; Vishniac 1983). The minimal mass of such satellites 
was estimated in Barkana, Haiman & Ostriker (2001). 

This means that the investigation of absorbers observed 
at high redshifts should be supplemented by the study of 
properties of dwarf isolated galaxies and discrimination be- 
tween such galaxies and dwarf satellites of more massive 
galaxies. It seems to be a perspective way to discriminate 
between models with the Jeans damping (58) and the damp- 
ing caused by the mass of WDM particles, (57), and, in par- 
ticular, between models with one and several types of DM 
particles. 
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